Binding Effects in Multivalent Gibbs-Donnan Equilibrium 
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The classical Gibbs-Donnan equilibrium describes excess osmotic pressure associated with confined 
colloidal charges embedded in an electrolyte solution. In this work, we extend this approach to 
describe the influence of multivalent ion binding on the equilibrium force acting on a charged rod 
translocating between two compartments, thereby mimicking ionic effects on force balance during 
in vitro DNA ejection from bacteriophage. The subtle interplay between Gibbs-Donnan equilibrium 
and adsorption equilibrium leads to a non-monotonic variation of the ejection force as multivalent 
salt concentration is increased, in qualitative agreement with experimental observations. 

PACS numbers: 82.60.Lf,82.35.Rs,82.39.Wj 



I. INTRODUCTION 

The classical Gibbs-Donnan equilibrium is usually invoked to explain various behaviors in ionic systems: swelling 
of polyelectrolyte gels, membrane potentials and osmotic shock experiments of cells or viruses pj. The physics behind 
these phenomena is associated with some colloidal charge (usually synthetic or bio-polyelectrolytes) confined within 
a subpart of the system, while small co- and counter-ions can diffuse freely in and out of this region 0] . Due to the 
electroneutrality inside the compartment, this confinement leads to some excess of counterions relative to the outer 
compartment, and therefore an additional osmotic pressure is set up inside. This effect is nowadays well-known and 
characterized experimentally on model systems Q. 

From a theoretical point of view, the Gibbs-Donnan approach is a simple and versatile way of incorporating 
electrostatic effects up to leading order in inhomogeneous systems, through the use of ionic chemical potential balance 
between two homogeneous phases and the electroneutrality condition. In systems either where colloidal charge is 
large or where the buffer solution contains multivalent ions, the Gibbs-Donnan theory has to be extended to take into 
account the electrostatic binding of counterions onto the colloidal charge [J] . Motivated by recent in vitro experiments 
measuring multivalent salt influence on the force ejecting DNA from bacteriophage Lambda H, we propose in this 
Letter such an extension in order to include multivalent ion binding effects in the classical Gibbs-Donnan approach. 
Indeed, we expect the interplay between Z-valent ion partitioning between inside and outside the viral capsid (pure 
Donnan effect), and the binding statistics onto DNA inside and outside (pure Langmuir-like adsorption statistics), 
to give raise to non-trivial translocation behaviors. In a slightly different context for example, the balance of these 
two effects has been shown by Klein Wolterink et al. to produce non-monotonic variation of the radius of gyration 
of a strongly charged polyelectrolyte star with respect to salt concentration Q. As will be shown in this work, this 
simple generalization of Gibbs-Donnan approach is able to explain at least qualitatively the non-monotonicity of 
ejecting force as function of multivalent salt concentration. In particular, we find that the repulsive force decreases 
and then increases, upon increase in the added salt concentration. The minimum in the force is associated to the 
charge neutralization point of DNA, which correlates well with experimental results pj. 

II. MODEL 

The model system considered in this work is depicted in figure^ It consists of a cavity of volume Vq embedded into 
a larger volume V >> Vq. The whole volume V is filled with a (Z:f ) electrolyte at an adjustable concentration. The 
cavity is permeable to both species (Z-valent cations and monovalent anions) of the electrolyte. A translocation gate 
on the cavity allows a rigid rod bearing L uniformly distributed negative charges to move between inside and outside 
the cavity. For a given configuration, there are Li n and L — Li n negative charges respectively inside and outside the 
cavity. Multivalent cations are expected to bind onto the charged rod, due to strong electrostatic interactions. The 
free energy of this system is written as F to t — Fi ayers + Ff ree + F neutra i where the first term is associated with the 
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FIG. 1: Sketch of the model. The two-sided arrows indicate different equilibria taking place in the system 

Z-cation adsorbed layers on the rod, the second term is associated to remaining free ions in the solution, and finally 
the last term takes into account any non-ionic effects. It can include for example the effect of a neutral osmotic 
pressure difference between inside and outside the cavity due to neutral polymers that cannot enter the cavity, as 
in the experimental s etup used by Evilevitch et al. [a, M\ , or bending effects in the case where the rod has a finite 

bending modulus [13, El El- 

Following the spirit of the original Gibbs-Donnan approach, the specific features of electrostatic binding are ne- 
glected, and the binding of Z-valent cations onto the rod is approximated by Langmuir-like adsorption isotherms . 
Each negative charge of the rod is therefore assumed to be a potential discrete binding site for a single Z-valent cation. 
The energy gain associated with a single binding event is given by —kTtA- This energy is of order ca ~ Zlg/b, where 
b is the typical size of rod unit and Ib the Bjerrum length. This neglects any correlation effects among adsorbed Z- 
valent ions. It can be checked a posteriori that including a certain level of correlations does not change the qualitative 
picture drawn by this simple model. Denoting by Naiu and NAout respectively the number of adsorbed ions on the rod 
inside and outside the cavity, the free energy of the adsorbed layers is Fi ayers = FA{NAin, Lin) + Fa{Naoui,L — Li n ), 
where the adsorption free energy of ions on a one dimensional lattice reads 

— = N A In — + {L- N A ) mil — J — N A e A (1) 

The first two terms describe the mixing entropy of adsorbed ions, and the last one is the adsorption energy. Assuming 
a unique value b 3 for the molar volume of each species (Z-cations, anions, and negative charges on the rod) for the sake 
of simplicity, the free energy of remaining free ions in the solution is given by the sum of translational free energies 
Ft(N, V) = kTN(\n(Nb 3 /V) — 1) of Z-cations and anions, respectively, inside and outside the cavity: 



7 free = <F T (N Zin - N Ain ,V - 6" Un) + F T (Z N Zm - L in ,V - b 3 Li, 
+ ^F T (N Z - N Zm - NAout, V-V - b 3 (L - L m )) 
+F T (Z(N Z - N Zm ) - (L — L m ), V-V - b 3 (L - L m ))\ (2) 



Note that the number of Z-cations inside the cavity is N Z i n , including both free and adsorbed ions, while the total 
number of Z-cations in the solution (both inside and outside the cavity, free and adsorbed) is N z . Following Donnan, 
electroneutrality conditions have been taken into account inside and outside the cavity, so that the numbers of inner 
and outer anions are AT_ in = ZN Zin - L m and iV_ ont = Z(N Z - N Zin ) - (L - L ln ). 

At equilibrium, the free energy of the system is minimum with respect to Naiu, NAout, ^z™ and The first three 
minimization conditions mean simply that the adsorbed layers are in equilibrium with the free Z-ions of the solution, 
both inside and outside the cavity (Langmuir balance), and that free Z-ions inside the cavity are in equilibrium with 
free Z-ions outside the cavity (Gibbs-Donnan balance). These equations are written in the limit V >> Vo, introducing 
species concentration n Zin = N Zin /V , n z = N z /V, (j) in = N A in/L in , (j> out = N Ao ut/(L - L in ), l in = L m /V Q , I = L/V 

i = (n Zin - (t)inlin)b 3 e eA (3) 
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FIG. 2: Ionic force acting on the rod vs rod length inside Li„, for different multivalent salt concentration. Parameters: 
L = 100, Vo = W^b^jV — l(fb 3 ,b — 2nm, Z — 4, eA = 2.5. Left panel: Total force; the large arrow and numbers highlight 
the variation of the force-length plots as multivalent salt concentration is increased. Note that the force decreases at first 
(1 — > 2 — > 3) upon increase in nz, and then increases (4 — > 5 — > 6). Upper right panel: Gibbs-Donnan contribution to the ionic 
force. Lower right panel: Langmuir contribution to the ionic force. 

— = (nz-J>outl)b*e^ (4) 



{nzin - 4>inkn) (Zn zin - l in ) = (n z - 4> out l) {Zn z - I) (5) 

For a given value of /,„ , these equations are solved self-consistently and their solutions allow calculation of properties of 
interest for the translocation problem, like the ionic contribution to the force acting on the rod fi on ic = d ( F,ay ° r ° +Ffree * > . 
Under equilibrium conditions, the value of U n is set by minimization of the free energy with respect to Li n . Introducing 
the force contribution associated with non-ionic features f neutral = dF j l £ ut, b al , this equation is simply the force balance 
on the rod fneutrai + fionic = 0. In this Letter, we focus on the behavior of the ionic force for different solution 
conditions at fixed Zj n , rather than determining the equilibrium partition of the rod inside and outside the capsid. 
Indeed, the value of U n is dependent on the choice of neutral force contribution j "neutral, which is beyond the scope of 
this work p|. 



III. RESULTS AND DISCUSSION 



In the absence of any analytical solutions for the equilibrium equations EE1 we solved them numerically for a given 
representative set of parameters. As it is seen in figure|5J the ionic force is always positive and increasing with respect 
to the rod length inside the cavity. This has to be contrasted with the net attractive force found in the work of 
Zandi et al. where the binding particles are only present inside the cavity and are not allowed to bind on the 
outer part of the rod. Now, increasing multivalent salt concentration leads to a non-monotonic variation of the force 
at fixed inner length in the present model, which is the main result of this work: the force is first decreasing at low 
salt concentration, and then increasing at higher salt concentration. The multivalent salt threshold corresponds to 
the neutralization point of the rod <pi n ,out ~ as it is seen on figure |31 This behavior is interpreted below by the 
interplay between Gibbs-Donnan and Langmuir contribution to the ionic force. 

The ionic force is rewritten as the sum of three main contributions fi on ic = Jgd + II + /aii- The three terms are 
respectively associated with the Gibbs-Donnan force contribution, the Langmuir force contribution and the osmotic 
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The first contribution Eq. [S] is associated with the free ions partitioning between inside and outside the cavity. In 
the absence of ion binding on the rod, there is a depletion of anions and an excess of Z-cations inside the cavity, due 
to the negative charge of the rod. Within classical Gibbs-Donnan approach, the net force is repulsive. The effect of 
ion binding on the rod at low multivalent ion concentration is mainly to reduce the net negative charge of the rod, 
without affecting the sign of the force. Above the neutralization threshold concentration n* z , for which (f> out ~ 1/Z, 
the net charge of the rod is positive, due to the binding of Z-cations. In this case, there is a depletion of free Z- 
cations and an excess of anions inside the cavity. However, in contrast to the previous case, there is still an overall 
excess of Z-cations inside (fig. EJ, a growing number of them being involved in the binding on the rod leading to its 
overcharging. The important proportion of adsorbed Z cations, not contributing to translational entropy inside the 
cavity, is now favorable to the rod being inside the cavity. As a consequence, the Gibbs-Donnan force is attractive as 
shown in figure This interpretation does not take into account the free energy of adsorbed layers, however, which 
compensates this negative Gibbs-Donnan force to give a net repulsive force. 

Indeed the second contribution, Eq. comes from the mixing entropy of adsorbed ions along the rod. The 
sign of the force is mainly given by the relative value of inner and outer adsorption degrees 4>i n and 4> ut- At low 
multivalent concentration, Gibbs-Donnan equilibrium favors free Z-cations excess inside, thereby increasing inner 
adsorption relative to the outside environment. Therefore the Langmuir force tends to pull the rod inside. This is 
the quasistatic effect seen in simulations of reference |8( for particular conditions, though with a larger amplitude, 
due to the asymetrical binding on the rod. Above the neutralization threshold, the situation is reversed, and outer 
adsorption is favored 4>out > 4>in ( n g- El) ■ As a consequence, Langmuir force is now repulsive, favoring rod ejection 
from the cavity. 

The last contribution Eq. [SI in the ionic force is the osmotic force, first introduced in [T^. and it is proportional to 
the osmotic pressure difference between inside and outside the cavity. It arises from the fact that, in order to insert 
the rod which has a finite volume inside the cavity, one has to perform a pressure work. Within the conditions chosen 
in this work, this osmotic force is one order of magnitude smaller than Gibbs-Donnan and Langmuir forces, and is 
therefore disregarded in the present analysis. 

For the chosen values of parameters (rod length, adsorption energy, etc.), the magnitude of the ionic force is of 
order of a few tenths of a piconewton, which is quite small at the molecular level. This comes mainly from the fact that 
the net entropic forces discussed previously arise from ionic concentration gradients between two compartments and 
two parts of the rod, unlike asymetrical situations where binding is limited to the inner part of the rod for example. 
Increasing the length and adsorption energy at fixed Vb produces larger ionic force amplitudes. Added to the neutral 
force contribution, the ionic force can still lead to substantial changing of behaviors, as is discussed below. 



As already mentioned in the introduction of this Letter, the present model was originally designed in order to 
evaluate the leading order of multivalent ionic effects on in vitro DNA ejection from bacteriophage 0. Within these 
experiments, DNA ejection is triggered using solubilized specific receptor in a solution containing neutral polymers 
producing an osmotic force that is able to balance the DNA ejection force. Changing ionic conditions of the solution 
at constant polymer concentration allows us to address the influence of ionic environment on the ejection force. The 
present model is of course a crude simplification of the real system. However it is able to predict qualitatively the 
non-monotonicity of ejection force, as it is observed in the relevant experiments. Within this model, the minimum in 
the force is associated to the neutralization of charged rod <fr ~ 1/Z. Using Eq. 21 the Z-valent threshold concentration 
scales like nz ~ e~ £A . This implies that this threshold concentration decreases with the valency of counterions, a fact 
that is observed experimentally for three different valencies: Na (Z = 1), Mg (Z = 2), spermine (Z = 4). Quantitative 
comparison of predicted and measured threshold is however difficult since the experimental buffer contains more than 
one multivalent salt. Our model can easily be extended to include more salt types and their binding competition with 
different energies CA i (i — 1,2...). This modified model leads qualitatively to the same results as the one derived for 
one multivalent salt (data not shown). The quantitative difference between the two models is the precise location of 



IV. APPLICATIONS AND LIMITATIONS 
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FIG. 3: Adsorbed and free Z-cation concentrations. The value of parameters are the same as in figure Left panel: Inner vs 
outer adsorbed degree of Z-cations for concentrations nzb 3 = 0.01,0.02,0.03,0.04,0.05,0.06 (cj>in,out increase with nz)-Upper 
right panel: Ratio of free Z-cation concentration inside/outside vs rod length inside. Symbols are defined in figure [5] Lower- 
right panel: Ratio of total Z-cation concentration inside/outside vs rod length inside. Symbols are defined in figure|21 

ionic force minimum or neutralization threshold. Since the main effect is already observed within the present simple 
one-salt model, we prefer to stick to this model, for the sake of clarity. Note that Gibbs-Donnan balance has already 
been used by Odijk and co-workers in order to describe DNA packing in a model-bacteriophage |14| . But the presence 
of multivalent ions as well as their influence on DNA ejection has been disregarded. 

The main effect that is not included explicitly within this model is the presence of correlations between bound ions, 
changing the effective adsorption energy |15|. However, we don't expect this effect to change drastically the results 
derived within our model: indeed, the non-monotonicity of the ionic force is related within our model to the presence 
of neutralization of the rod and possible overcharging, and not to its precise description. We checked with non-linear 
adsorption energy that prohibits the rod from being unrealistically too overcharged, that the qualitative behavior is 
unchanged. In particular, we observed that including only correlations and neglecting Gibbs-Donnan balance does 
not lead to non-monotonicity of ionic force, so that the experimentally observed minimum of the force can not be 
explained solely by invoking correlation effects. Therefore both basic ingredients of our model, Gibbs-Donnan balance 
and adsorption statistics, are necessary to observe the non-monotonicity of the force. Notice that for similar reasons 
we don't expect correlations of free ions, like the leading order Debye-Huckel correlations [l|, to change qualitatively 
the results obtained in this Letter. 

In order to interpret the results of DNA ejection experiments within this model, an additional implicit assumption 
has been made: the ionic effects can be decoupled from other contributions, termed "neutral" within this work. This 
assumption is valid within the simple geometry shown in figure ^ I n a more realistic model of viral bacteriophage, 
where DNA is arranged in a spool-like fashion within the viral cavity, the neutral contribution would include bending 
energy as well as direct electrostatic repulsion between neighbouring DNA turns in the spool 0, HTl Il2| . Up to now 
except in aforementioned works of Odijk and co-workers [l4(, the inclusion of electrostatic effects in such models has 
been based on phenomenological expressions describing osmotic pressure of hexagonal phases of DNA as measured in 
the experiments of [l6l | for exam ple. Since the simple extended Gibbs-Donnan model presented here is qualitatively 
consistent with these expressions 17f, we expect our model to describe correctly the leading order of multivalent ionic 
effects. 

The range of application of the simple model proposed in this Letter is larger than the ejection experiment just 
mentioned. Indeed it is relevant to the description of translocation properties of polyelectrolytes|jj, whenever the 
separating interface is permeable to mobile free ions. We showed in this paper that the force balance acting on 
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the rod can be displaced by changing ionic conditions. Within a biological context, the applicability of the present 
model for the translocation of bio-polyclectrolytes between different cell compartments, for example, is limited by the 
presence of ionic channels on the interface (cell compartment membrane) that actively regulate the ionic gradients 
accross the interface. Therefore the pure Gibbs-Donnan balance is no longer obtained, and new models have to be 
developed. 
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